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The cosmic microwave background (CMB) radiation provides a remarkable window onto 
the early universe, revealing its composition and structure. In these lectures we review and 
discuss the physics underlying the main features of the CMB. 


I. INTRODUCTION 

The cosmic microwave background (CMB) radiation provides a remarkable window 
onto the early universe, revealing its composition and structure. It is a relic, thermal 
radiation from a hot dense phase in the early evolution of our Universe which has 
now been cooled by the cosmic expansion to just three degrees above absolute zero. 
Its existence had been predicted in the 1940s by Alpher and Gamow mm and its 
discovery by Arno Penzias and Robert Wilson at Bell Labs in New Jersey, announced 
in 1965 [3] was convincing evidence for most astronomers that the cosmos we see today 
emerged from a Hot Big Bang more than 10 billion years ago. 

Since its discovery, many experiments have been performed to observe the CMB 
radiation at different frequencies, directions and polarisations, mostly with ground- 
and balloon-based detectors. These have established the remarkable uniformity of the 


E-mail: david.wands@port.ac.uk 
E-mail: oliver.piattella@pq.cnpq.br 
E-mail: casarini.astro@gmail.com 




2 D. Wands, O. F. Piattella and L. Casarini: Physics of the Cosmic Background Radiation 


CMB radiation, at a temperature of 2.7 Kelvin in all directions, with a small ±3.3 mK 
dipole due to the Doppler shift from our local motion (at 1 million kilometres per hour) 
with respect to this cosmic background. 

However, the study of the CMB has been transformed over the last twenty 
years by three pivotal satellite experiments. The first of these was the Cosmic 
Background Explorer (CoBE), launched by NASA in 1990 [H [3]. It confirmed the 
black body spectrum with an astonishing precision, with deviations less than 50 
parts per million [3]. And in 1992 CoBE reported the detection of statistically 
significant temperature anisotropies in the CMB, at the level of ±30 /rK on 10 degree 
scales [1]. COBE was succeeded by the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite, launched by NASA in 2001, which produced full sky maps in five 
frequencies (from 23 to 94 GHz) mapping the temperature anisotropies to sub-degree 
scales and determining the CMB polarisation on large angular scales for the first time. 
The Planck satellite, launched by ESA in 2009, sets the current state of the art with 
nine separate frequency channels, measuring temperature fluctuations to a millionth 
of a degree at an angular resolution down to 5 arc-minutes. Planck intermediate data 
was released in 2013 [7]0 

These lectures draw upon the excellent reviews of CMB physics by Hu and Dodelson 
IMI], Komatsu |12) and Crittenden |13) . We also refer the reader to comprehensive 
reviews on cosmological perturbations by Mukhanov et al. m and Malik and Wands 
m- Useful textbooks are those of Peebles m, Dodelson m, Mukhanov |18j 
and Weinberg |19j . Throughout this chapter we will use natural units such that 
h = ks = c = 1. 

II. BACKGROUND COSMOLOGY AND THE HOT BIG BANG MODEL 

We start by recalling the mathematical framework describing the expansion of the 
universe and the Hot Big Bang. Much of modern cosmology is based on general 
relativity and the framework of Friedmann, Lemaitre, Robertson and Walker in the 


^ After these lectures were given, full-mission data was released in 2015 [S], with Hnal polarisation 
data still to come. 
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1920’s and 30’s |20H22j . and Hubble’s discovery of the expansion of the universe 
We can “slice” four-dimensional spacetime into expanding three-dimensional space at 
each cosmic time, t, with a uniform matter density and spatial curvature. Requiring 
spatial homogeneity and isotropy at each cosmic time is known as the cosmological 
principle, which picks out the following space-time metric: 

dr‘^ 


ds^ = —df^ + a{t) 


+ r^dn^ 


( 1 ) 


1 — Kr^ 

where a{t) is the scale factor and k is the curvature of the maximally symmetric spatial 
slices, and we chose spherical coordinates with infinitesimal solid angle We will 
sometimes find it convenient to use conformal time, i], where dt = adr] and the line 
element takes the form 

dr'^ 


ds'^ = m(ry) 


-dr] + 


+ 


( 2 ) 


1 — 

The Hubble expansion rate is defined as H = d/a, where a dot denotes a derivative 
with respect to cosmic time t. The present value of H is called the Hubble constant 
and denoted as Hq. The value of Hq is often given in the form 


Hq = 100 h km s ^ Mpc 




( 3 ) 


Using the energy constraint, from Einstein’s equations of general relativity, one gets 
the Friedmann equation for the Hubble expansion 


= 


SttG 


A 

~P + ^ ~ 


3 " ' 3 a" ’ 

where we introduce the cosmological constant, A, and p, the energy density. The latter 
includes electrons, baryons (protons, neutrons and atomic nuclei), radiation (photons 
and neutrinos) and dark matter (non-baryonic massive particles, non-relativistic by 
the present day). 

Dividing through by , eq. Q can be cast in the following dimensionless form: 


1 — -j- “t“ . 

where we define the relative contributions to the Hubble expansion 

SirGp „ A 


n = 




Da = 


3id2 ’ 


= _ — 

“ a2dd2 


( 5 ) 


( 6 ) 
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In order to get a closed system of equations we must determine the evolution of the 
density, p in eq. as a function of the scale factor. For this we can use the continuity 
(energy conservation) equation 

p = -^H{p + P), (7) 


plus an equation of state for the pressure, P{p). We will be interested in three 
important cases: 

• n = rir = 1, radiation domination: 

Pr = -Pr => Pr OH O OC OC 7/ . (8) 

o 

• = Clm = 1, matter domination (Einstein-de Sitter): 

Pm = 0 Pm OC a~^ a OC OC if. (9) 


• Oa = 1, a domination (de Sitter): 

a OC OC {poo - p)~^ 


( 10 ) 


The CMB consists of photons which survive from an early, radiation-dominated, 
Hot Big Bang and have a small density with respect to non-relativistic matter today. 
Nonetheless the CMB holds a rich store of information about the history of our 
Universe, as we shall see. For example, recent observations of the CMB by Planck 
|24j can be used to infer values for the above cosmological parameters at the present- 
day: 

h = 0.674±0.014 , flo = 0.314±0.020 , IIao = 0.686±0.020 , = -0.04±0.05 . 

( 11 ) 

The data are consistent with a flat universe, k = 0, which will be our working 
hypothesis hereafter. We see that the expansion today is dominated by a cosmological 
constant (or some form of matter which acts very much like a cosmological constant) 
but in the recent past it was dominated by non-relativistic matter, and before that by 
radiation. 
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A. Black-body spectrum 


The CMB is observed to have a black-body spectrum characteristic of a thermal 
equilibrium distribution, consistent with the hypothesis that our Universe emerged 
from a hot, dense Big Bang. 

Photons follow a null trajectory in the FLRW metric © such that 



( 12 ) 


where h* is a unit 3-vector, gijh^h^ = 1. The 3-momentum of a photon is p* = ph*, 
where p is the wavenumber (remembering that we are using units such that h = 1 and 
c = 1, so that p also describes the energy of a massless photon). 

CMB photons have an isotropic Bose-Einstein distribution function with temper¬ 
ature T 


fip) = 


1 


(13) 


exp(p/T) — 1 

Given this isotropic distribution, we can compute the number density of CMB photons 


riy = 2 




2.4. 


vr^ 


(14) 


where the photons have 2 independent polarisation states and Anp^dp is the volume of 
an infinitesimal shell in three-dimensional momentum-space. Their energy density is 

Airp^dp 




( 15 ) 


However the CMB photons are no longer in equilibrium with the matter we see 
in the universe today. The photons are free to propagate through the universe after 
electrons and baryons have recombined into neutral atoms, so the black-body spectrum 
must be propagated to the present day from the early universe. Freely propagating 
photons follow the geodesic equation in curved space-time 


dP^^ 

hx 


+ = 0 , 


(16) 


where T^o- is the Christoffel symbol. We define the photon 4-momentum as P^ = 
dx^/d\, where A is an affine parameter, and the modulus-squared of the 3-momentum 
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is = pijP^P^ where gij is the spatial part of FLRW metric Q. From the geodesic 
equation in the conformal FLRW metric Q we obtain 

1 dp 1 da , 

p dr] a dr] 

Integrating this up to the present we obtain the cosmological redshift of the photon 
momentum, defined as 


, , _ P ao 

1 + z = — = — 
po a 


(18) 


We can interpret this simply as the expansion of the universe stretching the wavelength 
of a photon, reducing (redshifting) its energy and momentum. 


Note that the form of the Bose-Einstein distribution (13) is preserved 

1 


f{p) = 


1 


exp(p/r) - 1 exp(po/Eo) - 1 ’ 
where the temperature is also redshifted with the expansion 

T 

l + z = — . 

To 


(19) 


( 20 ) 


Thus we see that the energy density (15) of the photons decreases as the universe 
expands 

.-4 


P-y 


oc a 


( 21 ) 


Although photon density is small is in the universe today, it dominated the hot, dense, 
early universe. 


B. Hot Big Bang 

At sufficiently high temperatures we expect all particles to be relativistic. If these 
particles interact and efficiently redistribute energy they will share the same thermal 
equilibrium temperature. To be relativistic we require T ^ m, i.e., the thermal energy 
is much larger than the rest mass of a given particle species. At this stage of the 
primordial universe we can write the energy density using the same form given in 


eq. (15) for all the relativistic species: 








The Cosmic Microwave Background Radiation 


7 


where is the sum of the effective number of degrees of freedom. Each bosonic 
species in thermal equilibrium contributes one per spin state (e.g., photons contribute 
+2, corresponding to two polarisations), whereas each fermion contributes 7/8 per spin 
state, due to the different statistics 0 

In a radiation-dominated universe Q the time dependence of the scale factor is 
given by a oc t^/'^ and thus from eq. Q we have 

3/^2 3 


P = 


SttG 327rG'f2 


so that from (22) time and temperature are related by 

1 


t = 


30 


327rGg,s7T^T^ ' 

Thus we have the simple, approximate temperature-time relation 

t _ 1 /iMeVy 

1 sec \ T J ' 


(23) 


(24) 


(25) 


C. Spectral distortions 


The black-body shape of the CMB spectrum is maintained at early times because of 
the high interaction rate of photons with the other particles of the primordial plasma. 
We can identify two principal scattering processes which contribute to maintaining an 
isotropic, equilibrium distribution: 

• Compton scattering: scattering of photons and relativistic electrons, redistrib¬ 
uting energy and momentum, conserving photon number 

e“ -|- 7 -H- e“ -|- 7 . 

At low energies this reduces to Thomson scattering, i.e., elastic scattering of 
photons off non-relativistic electrons, exchanging momentum, but conserving 
photon energy and number. 

^ If a species decouples from this thermal bath, but remains relativistic, it can contribute with a 
different temperature in the above equation. This is what happens for neutrinos. They decouple 
relativistically from the primordial soup, at T « 1 MeV and their temperature today is expected to 
be (4/11)^/® times that of the photons because photons are heated by annihilation. 
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• Double (radiative) Compton scattering: scattering of photons and relativistic 


electrons, redistributing energy and momentum, and changing photon number 


e d-y-H-e +7 + 7 . 


Many processes in the early universe before the time of recombination could 
potentially lead to measurable distortions in the CMB spectrum, which might be 
measured with future missions. Particle annihilation or decay would heat the 
primordial plasma, and hence the photons, or even the evaporation of primordial black 
holes in the relevant mass range. Even the damping of small scale density variations 
in the primordial plasma due to photon diffusion can lead to deviations from an exact 
black-body spectrum. For more detail about CMB spectral distortions and what might 
cause them, see [25] . 

Efficient Compton and double Compton scattering maintains a full thermal 
equilibrium spectrum above a redshift [ 10 ] 



(26) 


where determines the density of baryons and hence (in an electrically neutral 


universe) electrons. 

Below this redshift Compton scattering can still redistribute energy and momentum 
between photons and electrons, but double Compton scattering becomes inefficient. 
In the absence of double Compton scattering, interactions cannot create or remove 
photons from the plasma. Compton scattering still maintains a statistical equilibrium 
above redshift [TO] 



(27) 


Thus if additional energy is dumped into the primordial plasma below redshift zth the 
CMB photons acquire a statistical equilibrium distribution 


f{p) 


1 


exp[(p - n)/T] - 1 ■ 


(28) 
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with non-zero chemical potential /r. This is known as a //-distortion in the CMB 
spectrum. Limits from the COBE satellite give an upper limit on the size of such a 
distortion [6]: 

^ < 9 X lO-*^ at 95% CL. (29) 


Below the redshift Compton scattering off relativistic electrons becomes ineffi¬ 
cient. High-energy electrons along the line of sight can still transfer energy to low- 
frequency photons via inverse Compton scattering, without reaching statistical equi¬ 
librium. This leads to a characteristic “y-distortion” where low energy photons are 
boosted to higher frequencies, leading to a deficit in the CMB intensity at low frequen¬ 
cies in the Rayleigh-Jeans region, equivalent to a temperature deficit 


AT 


p<.T 


-2y. 


(30) 


and an enhancement at high frequencies. The Compton //-parameter is defined as the 
line-of-sight integral of the electron pressure 


y = 


— —neC^Tdl ■ 
me 


(31) 


where rig is the density of free electrons and or is the Thomson scattering cross-section, 


see Eq. (33) below. Constraints from COBE/EIRAS give the upper limit [6] 


\y\ < 1.5 X 10“^ at 95% CL. 


(32) 


These constraints still rely on COBE observations, more than twenty years ago. 

An important source of //-distortions seen in specific directions in the CMB is 
the Sunyaev-Zeldovich effect [26], from hot cluster gas along the line of sight after 
recombination. The Planck satellite has now compiled a catalogue of 439 clusters 
detected in the Planck data via their SZ signal m with many more being detected 
by ground-based experiments such as the Atacama Cosmology Telescope [28] and the 
South Pole Telescope |29|. 
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D. Tight-coupling and sudden recombination 


At low energies (much smaller than the electron rest mass) electrons and photons 
interact via Thomson scattering, whose cross-section is 


(Tt = 


3m2 


= 6.65 X 10“^® m^ 


(33) 


The corresponding mean-free-path for photons associated with Thomson scattering is 
given by 


^mfp — • (34) 

rieO-T 

Around z ~ 1100 the mean-free-path is approximately 2.5 kpc, corresponding to a 
comoving scale of order 2.5 Mpc at present |10) . On scales much larger than the 
mean-free-path, A ;§> Amfp, the photons are tightly coupled to the electrons, while 
electrons are tightly coupled to protons through the Coulomb interaction. In this 
regime, photons, electrons and protons can be treated as a single fluid with common 
3-velocity, and isotropic pressure. 

The mean-free-path is time-dependent because the free-electron density, rie, is time- 
dependent. As the Universe cools down the capture of electrons by protons becomes 
efficient. As the wavelengths of photons are redshifted by the cosmic expansion, fewer 
photons have sufficient energy (the ionisation energy, 13.6 eV) required to break the 
binding energy of an electron in a neutral hydrogen atom. Therefore, the density of 
free electrons, Ue, rapidly drops around z « 1100, leading to a rapid increase in the 
Thomson mean-free-path beyond the Hubble radius. 

This process is called decoupling, because photons no longer interact with 
electrons. It is also called recombination because this is the epoch when protons 
and electrons recombine to form hydrogen atoms. Recombination and decoupling are 
practically simultaneous because the rapid drop in the density of free electrons due 
to recombination affects the Thomson scattering rate. By solving the corresponding 


^ The full cross-section describing the process e“ -I- 7 —>■ e“ -I- 7 is given by the Klein-Nishina formula 
|30 |. which displays not only the dependence on the photon energy but also on its polarization and 
the scattering angle. Since the energies involved in the recombination process are much smaller than 
the electron mass, we can safely use Thomson cross-section. 
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Boltzmann equation we see that recombination and decoupling occur at redshift |10] 

/ r\ 1,2 \ 0.0105 / 7 2 \ — 0.028 

l + .. = 108!>(5=t!l (iW!) , (35) 

V 0.14 ) \o.024j ’ 

Note that this is some time after (but not long after) matter-radiation equality, 


1 + Zeq = 3.4 X lO'" 


0.14 


(36) 


Another way to define when recombination/decoupling takes place is via the 
Thomson optical depth 

rvo 


r = / UeUrdt , 
Jr] 


(37) 


which represents the integrated scattering rate from a conformal time r] until today 
rjQ, i.e., the average number of scattering events between these two times. The spatial 
hyper-surface of constant rj = where 77 * is the conformal time corresponding to r = 
1, is called the last-scattering surface. Of course, recombination is not an instantaneous 
phenomenon, but it occurs sufficiently rapidly that a useful approximation on comoving 
scales greater than about 2.5 Mpc is the so-called sudden recombination, as if it really 
happened at a single instant, 7]^,. 

III. CMB ANISOTROPIES 


Anisotropies observed in the CMB radiation are caused by inhomogeneities in the 
cosmological spacetime and matter distribution. Fortunately, these inhomogeneities 
are small (about 1 part in 10^) with respect to the background homogenous energy 
density, thereby allowing us to use perturbation theory to model their behaviour. In 
the following we shall consider a linearly perturbed distribution. 

We do not measure the plasma density directly, but rather anisotropies, in the CMB 
photon distribution function, /—)■/ + Sf. At first order these can be described by a 


perturbation in the temperature of the Bose-Einstein distribution (13), where 


T{r], X, p) = T{r]) [1 -k ©(r?, x, p)] 


(38) 
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where p denotes the direction of the photon propagation. The temperature fluctuation 


in the plasma is related to the photon density contrast via Eq. (15) as 


_ 5T _15p^ _ 1 


(39) 


A. Spherical harmonics 


Since we observe CMB on the celestial sphere, it is useful to expand 0 in spherical 
harmonics 

00 ^ 


^ ^ a£m(?/,x)y£m(p) • 


(40) 


£=0 m=—£ 

since the spherical harmonics form a complete orthonormal basis on the sphere 


I dnnYim{n)Y;,^,{n) = 6ee5mm' . (41) 

The coefficients aim describe the temperature fluctuations at a given angular multipole 
1. An isotropic distribution has an angular power spectrum Cf. 




(42) 


In this case the correlation between the temperatures in two directions on the CMB 
sky depends only on the angular distance between the two directions and not on the 
orientation of the arc which joins them. 

For a fixed i, one has 2£ + 1 different be., 2£ + 1 independent estimates of 

the true C^. The “observed” corresponds to our best estimate of the true angular 
power spectrum: 




* ^obs 
^Im 5 


(43) 


i.e., it is an average over the observed multipole moments, m, at fixed i. We define the 
cosmic variance as the expected error in our determination of the true power spectrum 


ACi 

Cf 


cosmic variance 


'Ci-c: 

Cf 


obs' 


(44) 
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Calculating the expectation in the above equation, with the help of eq. (42), one obtains 


AC, 

Q 


cosmic variance 


21 + 1 


(45) 


Thus at small multipoles, £, corresponding to very large angular scales, the cosmic 
variance is significant and represents the minimal uncertainty in estimating the true 
angular power spectrum given that we have only one realisation of the CMB sky. 


B. Last-scattering sphere 


Since most photons are last scattered at r/*, we will be mostly interested in their 


distribution, 0 (t/, x,p) in Eq. (38), at evaluated at recombination, i.e., at initial 
time r] = r]^ and comoving displacement with respect to an observer at the origin, 
X* = — Zl*p, where the comoving distance to last-scattering = r/o — — r/o. Then 

we propagate this photon distribution until today using the free-streaming equations, 
i.e., the collision-less Boltzmann equation for photons. 

Adopting the sudden-recombination approximation, we assume that the photons 
are tightly coupled with an isotropic distribution up until last scattering. 


0*(p) = 0 (r/*,x*) . (46) 

The CMB temperature varies across our sky due to the variation in the photon 
temperature across the last-scattering surface. 

We can decompose this 3D CMB temperature field into Fourier modes 

0(r/,x) = I d^k0(r/,k)e**^- . (47) 

Linear modes with different comoving wavevectors, k, then evolve independently at 
first order. We assume that these perturbations are stochastic quantities drawn from 
some distribution, which usually is assumed to be Gaussian. 

The expectation value of each mode is zero and its variance is the power spectrum 


(0*(T/,ki)0(r/,k2)) = (27r)^(5^(ki + k2)P0(A:i,r7) . 


(48) 
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Note that Pq is function of the modulus of ki only, i.e., we assume statistical isotropy. 
The correlation function in real space is given by the Fourier transform of the power 
spectrum 

1 


^ 0 (r) = ( 0 ( 77 ,x) 0 (T/,x + r)) = 


d'^kP'^’^Peik) . 


(49) 


(27r)3 j 

Angle brackets denote the ensemble average. That is, one imagines different 
possible realizations of our universe. In theories such as inflation, where primordial 
fluctuations are quantum in their origin and then become effectively classical through 
an exponential phase of expansion, it is possible to predict the primordial form of the 
power spectrum. After that, it is evolved up until today using the classical equations 
of cosmological perturbation theory. Thanks to the ergodic theorem, we can swap the 
ensemble average into a position average, see Appendix D of [T^ . 

Since Pq depends only on the modulus k, we can perform the angular integration 


in (49) and find 


^e{r) = 


27r^ 


''°°dk n sm kr 

-k'^Peik)- 


'0 


k 


kr 


From the above result, we can identify the dimensionless power spectrum 

_ k^Peik) 

Pe{k) = . 


(50) 


(51) 


We can decompose the temperature field on the last-scattering surface into spherical 
harmonics using the plane-wave expansion 

OO I 


gik-r ^ 


1=0 rn=—l 


(52) 


where the spherical Bessel function ji{x) is defined in terms of the regular Bessel 
function Ji-\-i/ 2 {x) as j^{x) = {-k/ 2xY^‘^ Jij_ii 2 {x). Substituting this expansion into (47) 
and comparing with (40) evaluated at x* = — we obtain the spherical harmonic 
coefficients 




(53) 


and hence the angular power spectrum (42), by using eqs. (41) and (48), becomes: 

/ OO nu 

—VQ{k)jl{kD,) . 


(54) 
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The window function 


Wt{k) = A7rj,{kD,)^ 


(55) 


peaks about k = i/so one obtains approximately that 


Ve{^ho) , 


(56) 


by using « ryo and the result 



2l{l + l) ■ 


1 


(57) 


This is the origin of the ubiquitous prefactor l{l + 1) in CMB spectrum plots. In order 
to obtain the full result one should include contributions from the metric perturbations 
and the dipole at recombination and the ISW effect, which we present in the following 
section. 


IV. SACHS-WOLFE FORMULA 


In the previous section we discussed the basic quantities which describe the 
CMB temperature anisotropies at last-scattering, and in particular the angular power 
spectrum, C^. In this section we link these to the observed temperature fluctuations 
including the effect of inhomogeneities in the metric and the density distribution of 
the matter content in the universe. We will derive the Sachs-Wolfe formula. In order 
to do this, we present the essential elements of relativistic cosmological perturbation 
theory, focusing on first-order fluctuations. The pioneering work in this field is due to 
Lifshitz m but we also refer the reader to more recent reviews, such as m- 


A. Metric perturbations 


The starting point for discussing cosmological perturbations is the perturbed FRLW 
metric |15j 


ds'^ = { — (1 + 2A)drf‘ + 2V iBdx’^drj + [(1 + 2C)6ij + 2Vi'VjE] dx^'dx^'^ , (58) 
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where A, B, C, and E are scalar functions of the coordinates. In the above metric, 
we are considering only scalar perturbations, neglecting for now vector and tensor 
(gravitational wave) perturbations. Because of the tensorial nature of the metric, the 
above scalar functions change when changing the reference frame. It could happen 
that a reference frame exists in which A = B = C = D = 0. In this case then there 
are no metric perturbations, since we recover the original unperturbed FLRW metric. 
So, the fact of having four scalar functions of the coordinates in the above metric does 
not guarantee that we are actually dealing with cosmological perturbations, because 
the latter may be coordinate artifacts. This is the well-known gauge problem. 

In order to know if we are really dealing with cosmological perturbations, a useful 
tool is to construct combinations of the above scalars which remain invariant under 
first order coordinate changes. There are three combinations independent of the spatial 
threading: A, C and a = E' — B^ where the prime denotes differentiation with respect 
to the conformal time rj. There are then two combinations independent of time slicing, 
for example, the Bardeen potentials [laiisiES] 

^ = A-na-a', ^ = C-fia . (59) 


In the above definition % = a'/a, is the conformal Hubble parameter, i.e. defined with 
respect to the conformal time. 

A particularly useful gauge is the conformal Newtonian gauge, where the metric 


becomes diagonal since the choice is B = E = Q. The Bardeen potentials (59) can be 
identified with the metric perturbations A and C in this conformal Newtonian gauge 
(where a = 0). The perturbed metric thus takes the form [TO] 


ds^ = { — (1 + + (1 + 2^)6ijdx''dx^^ . 


(60) 


It can be shown by writing down explicitly the Einstein equations that their spatial 
traceless part depends on <I> -|- \I'. For example, the quadruple moment of the matter 
distribution acts as source of the spatial traceless part of the Einstein equation. In the 
tight coupling limit, there is no anisotropic stress because the high interaction rate of 
photons due to Thomson scattering establishes an isotropic distribution, which implies 
that <1> -I- T = 0. 
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One can construct other gauge-invariant variables, e.g., involving matter quantities, 
such as the density contrast and velocity potential in the conformal Newtonian gauge 

5p — p' a 


5 = 


P 


V = v + E' 


(61) 


or the curvature perturbation 


C^C-^5p, (62) 

P 

which can be identified with the metric perturbations C in the uniform-density gauge. 
This is a particularly useful variable on large scales since ^ is conserved for adiabatic 
perturbations on super-Hubble scales {k <C aH) [33]. For example, simple slow- 
roll inflation models typically produce an approximately scale-invariant dimensionless 
power spectrum, V^{k), on large scales at the start of the radiation dominated 
era. Thus we will typically set initial conditions in terms of C and/or isocurvature 
perturbations. 

Note that these different perturbation variables are not necessarily independent. 
For example we can express ^ in terms of the conformal Newtonian gauge quantities: 

C = ch-^5. (63) 

P 


B. Perturbed geodesics 


What is the form of a perturbed geodesics in the conformal Newtonian gauge (60)? 
By setting ds^ = 0 for a null trajectory, we find the coordinate velocity of a photon 




(64) 


where p* is a unit vector, = 1. Defining the 4-momentum as = dx^/dX and 

the modulus of the 3-momentum = pijP^P^, the perturbed geodesic equation (16) 
can be written as follows: 


1 dp /1 da \ ST 

p dr] \a dp dp J ^ dx^ 


(65) 


The term in parenthesis is the usual Hubble redshift corrected by the metric 
perturbation, which makes the expansion not homogeneous and isotropic, as it was 
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in the background. The last term represents the gravitational blueshift or redshift 
experienced by a photon falling into or climbing out of a potential well. Introducing 
the total time derivative along the photon path, i.e. 


dr] dr] ^ dx^ ’ 


the geodesic equation (65) becomes 


1 dp 1 da 

pdr] a dr] 


d'k d 


( 66 ) 


(67) 


This can be formally integrated along the photon trajectory from recombination, r/*, 
until today, r]Q, 


In = _ In - ^0 + (^' - . (68) 

Splitting the momentum in a background part plus perturbation, i.e. p —>• p + 6p, one 
obtains 


6p 


6p 


rvo 


^ = — I + - ^'o + 

P Jo \P 


(^' - ^')dr] 


(69) 


'V* 


This relative perturbation in the photon momentum causes a relative temperature 
fluctuation in the CMB, 0 = 6p/p. So, one sees that at recombination photons get 
a redshift escaping from over-densities on the last-scattering surface with a negative 
gravitational potential 'k*. This is part of the Sachs-Wolfe effect [Mj- The integral term 


in (69), named the integrated Sachs-Wolfe effect, accounts for the time-dependence of 
the potentials along the line of sight from recombination until today. 

The observed temperature fluctuation includes a Doppler shift due to the relative 
peculiar velocity (in addition to the expansion) between the last-scattering surface and 
the observer 


0=^ + n.V. 

P 


(70) 


The full Sachs-Wolfe formula is thus 


I rvo 

©obs = -k T* - h • V* -k / (^'' - ^')dr] - To -k fi • Vobs , 


( 71 ) 
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where we have identified the relative momentum perturbation for photons on the last 
scattering surface with the radiation density contrast 5-y = ^b’pj'p. The first three 
terms represent the intrinsic Sachs-Wolfe effect (on the last-scattering surface) and the 
fourth the integrated one we already mentioned, iho is the gravitational potential at 
the observer today and gives an undetectable correction to the monopole (that is, the 
solid-angle-averaged temperature). The last term is a dipole anisotropy induced by 
the observer’s velocity. 


C. Adiabatic and isocurvature perturbations 


In order to evaluate the relative contribution of different terms in the Sachs-Wolfe 
formula ( |71[ ) we need to determine the evolution of linear perturbations, in particular 
at the time of last scattering. 

The behaviour of the scalar perturbations previously introduced is given by the 
Einstein evolution equations (coming from the spatial part of the Einstein equations 
written in an arbitrary gauge): 


C + 2UC' 


- (th! A = 
a + lUa -C- A = 


-AttGo^ (^5P + , 

(72) 

SyrGa^n , 

(73) 


subject to the Einstein energy-momentum constraints (the time-time and time-space 
components of the Einstein equations): 

m{-C + UA) + V2(C - Ua) = -^T^Ga^dp , (74) 

-C+ 'HA = -^T:Ga^{p + P){v + B) . (75) 


Einally, not independent from the above equations, we also have the energy and 
momentum conservation (continuity and Euler) equations: 


5p' + m{5p + 5P) + 3(p + P)G' + {p + P)V\v + E') = 0, 


(v + py + (1 - 3ci)niv + B) + (j) + 


p + p 


5P+-V^U 1=0. 
o 


(76) 

(77) 


Note that the fluid quantities above introduced {6p, 6P, etc.) refer to the total matter 
content, but if the components do not interact among themselves, these equations 
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can also be considered individually for each one of the components that make up the 
balance of the cosmic energy budget. 

Exploiting the gauge freedom, we may consider the continuity equation written in 


the uniform density gauge, i.e. for Sp = 0 where C = C* in Eq. (62): 


3?^<5Pnad + 3{p + P)C' + ip + P) W = 0 , 


(78) 


where we identify the pressure perturbation with the non-adiabatic pressure in this 
uniform-density gauge, (JEand = 6P— {P'/p')5p, and the conformal Newtonian velocity 


V was defined in Eq. (61). Rearranging (78) we have 

^ <^^nad _ . 


(79) 


For fluids with a barotropic equation of state, P = P{p), we automatically have zero 
non-adiabatic pressure, 5Pnad = 0. Thus, on large scales, where the contribution from 
the divergence of the conformal Newtonian velocity, V^E, can be neglected, we have 
C being conserved. 

The same argument can be applied to any non-interacting barotropic fluids 


Thus, generalising the definition of the curvature perturbation (62), we get conserved 
perturbations on large scales for radiation and matter, (^.y and Cm- These can be written 
in terms of conformal Newtonian gauge quantities as 

< 5 . 


C. = ^ + $, 

Cm = ^ + <h 


(80) 

(81) 


Initial conditions are set up at sufficiently early times and on very large scales. 
Considering only radiation and matter, we have the total curvature and entropy 
perturbations 

= V 7 C, + 4pU.. _ _ 3( _ j ^ 

Apy +3pm 

Adiabatic initial conditions are defined as: 


C = Cm = C7 = constant 


5’m = 0. 


(83) 
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And the isocurvature initial conditions are defined as 


C-y = 0 , Sm = 3Cm = constant . 


(84) 


A full treatment of the initial conditions requires neutrinos and baryons to be dealt 
with separately, giving rise to two more isocurvature density modes [36]. 

The above constants are in general dependent on the scale. This dependence 
is set during an inflationary era and is thought to come from primordial quantum 
fluctuations. A slow time-dependence of the evolution during inflation leads to a weak 
scale-dependence of the dimensionless power spectrum. 


dlnV^ 


(85) 


The intrinsic Sachs-Wolfe effect ( (7l| ) on large scales (neglecting the velocity V*) 
can then be written as 


^ _ 1 

~T ~ 


-|- 4'^, — C 7 T 24* 


( 86 ) 


Let’s now see how the evolution of the gravitational potential looks like, since we have 
seen its role in determining the Sachs-Wolfe effect. Consider a barotropic fluid with 
equation of state p = wp and constant w. Neglecting anisotropic stresses, so that 
<f> = — T, the evolution equation is computed from the spatial trace of the Einstein 


equation (72) and reads: 


T" -h 3(1 -k = 0 . 


(87) 


On large scales or during matter domination (to = 0), where one can neglect the spatial 
gradient, we obtain a constant solution 


To = - 


3{l + w) 
5-k3'«; 


( 88 ) 


where we related this constant to C which we have already shown to be conserved on 
large scales and for adiabatic perturbations. This sets the initial conditions for the 
scales which subsequently enter the horizon. 

The evolution of the gravitational potential behaves in a two very different ways 
from the radiation dominated phase compared to the matter dominated era. In 






22 D. Wands, O. F. Piattella and L. Casarini: Physics of the Cosmic Background Radiation 


radiation, w = 1/3, for a comoving wavenumber k, we find 


^kiv) 


C7 


6 

{kr]y 


cos 



6v/3 . 
{kr])' 


sm 



(89) 


For super horizon scales {kr] <C 1) the above solution tends to a constant: —)■ 

—2(^/3. But for sub horizon scales {kr] —)■ 00 ) the potential oscillates and decays, 
'hfc —)> 0. Thus the growth of matter inhomogeneities will also be suppressed in this 
regime. 

In the matter dominated era {w = 0) the gravitational potential is constant at all 
scales 


= -3Cm/5. 


(90) 


With this result, the intrinsic Sachs-Wolfe effect on large scales can be written as 


6T 

~¥ 


C 7 T 2 '!'^, 


C7 g Cm ) 


(91) 


where we used the matter-dominated solution for the gravitational potential because, 
as previously noted, recombination takes place in this regime, z* < Zgq- For adiabatic 
perturbations {Sm = 0 and Cm = C7) the contribution is: 


5T 

~¥ 


ad 


= --C-7 = -'k*, 

5 S 7 3 * ’ 


(92) 


whereas for isocurvature perturbations {(^ = 0, Sm = 3Cm): 


5T 

Y' 


= -zSm = 2T* . 
5 


(93) 


With these two formulas, the large-scale approximation (56) can be written as follows: 

l{l -|- 1 ) 


27r 


-Cl ~ ^'P({i/r]o) + i^'Ps{iho) + ^C(:s{i/'no) 


(94) 


i.e. a contribution coming from the adiabatic perturbations, another from the 
isocurvature perturbations. Vs, and their cross-correlation, C(^s (recall that the power 
spectrum is a quadratic function of the perturbation variables). 
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V. ACOUSTIC OSCILLATIONS 

One of the striking features of CMB power spectrum is the presence of acoustic 
oscillations, which originate from sound waves in the baryon-photon fluid at the time 
of last scattering. We present here the physics which lies behind this phenomenon, 
following the approach of Hu m- 

The coupled first-order energy and momentum conservation equations for radiation 
perturbations in conformal Newtonian gauge (neglecting for the moment the effect of 
baryons) are 

= (95) 

V; =VT . (96) 

Eliminating yields an oscillator equation: 

+ '^) - 3^' (4'^7 + 'k) = ('h - CD)" . (97) 

A. Matter era 


It is easy to solve this equation in the matter era, since T = —‘h = constant, 
therefore the right-hand-side vanishes and we get 

^ cos(/cs) , (98) 

for k ;§> keq, where keq is the comoving wave-number of a scale which crosses the 
Hubble horizon at matter-radiation equivalence, i.e. kgq = i^eqOeq- The sound horizon 
for the relativistic fluid is given by 


■s = ^ Csdrj' = -^Tj . (99) 

We shall see later that the presence of baryons affects the sound horizon since it 
modifies the speed of sound. 

At recombination, for adiabatic perturbations, the above solution can be written 


as 



1 

5 


C 7 cos(/cs*) . 


( 100 ) 
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We see in Fig. the characteristic oscillating behaviour of the CMB temperature 
fluctuations, known as acoustic oscillations. Note that the power spectrum corresponds 
to the square of the amplitude so that peaks in the power spectrum occur at maxima 
or minima of the amplitude. 



Figure 1. Evolution of the effective temperature, eq. (98), and its absolute value (corresponding 


to the power peaks, solid lines) PITU|. 


The first peak takes place when the argument of the cosine is equal to tt, i.e. the 
comoving acoustic scale, Xa, is defined as 

kAS* = TT =► . (101) 

V3 

Dividing by the comoving distance to recombination, we find the angular scale 
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which in the matter-dominated universe can be approximated as 

— ^ ^-1/2 ^ 2° . (103) 
VO 

This is the angular scale of the particle horizon at recombination. It spans only 2 
degrees in the sky, and yet we see a high degree of isotropy in the CMB sky on all 
angular scales. No causal process could lead to isotropy on scales separated by more 
than 2° in the classical hot big bang because of the lack of casual connection on any 
large scales. This is the well-known horizon problem in the big bang model which is 
most clearly seen in the CMB, and to which inflation offers a solution. 


B. Radiation driving 


In the radiation dominated era, the gravitational potentials cannot be assumed to 
be constant and they rapidly oscillate and decay as each scale enters the horizon, i.e. 
$ = —T—)-0asfc7y—>^oo. Thus we can neglect the source term in the oscillator 


equation (97) for A:r/ » 1 to obtain 


-(5.y -|- 'k j ~ -<5-y oc cos(fes) 


(104) 


Extrapolating forward to the matter era, we obtain the same oscillations for the 


radiation density as before, Eq. (100), but with larger amplitude the potentials are 


negligible. For scales that enter the horizon well before matter-radiation equality, 
k ^ keq, we find 


1 


-|- T I cos{ks^:) 


(105) 


Comparing with Eq. (100), we find that the amplitude of oscillations for k > keq can 


be five times larger for low matter density, as this delays matter-domination, which 
causes the gravitational potential to decay on sub-horizon scales during the radiation 
era, see Figure 
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Figure 2. Evolution, starting in the radiation era, of the gravitational potential, 'h (red), the 
effective temperature, (1/4)5^ +(blue), and the anisotropic pressure due to photon diffusion, 
TT^ (green) [TU] , 


C. Baryon loading 

What happens if we include the effect of baryons which are strongly coupled to 
the electrons, and hence photons, before recombination? Intuitively, since baryons are 
massive particles, they would slow down the oscillations. 

We define the coupled baryon-photon velocity, Vfyy, using 

(p-y + P^)\.y + (ph + Pi|)Vb = (1 + R){p^ + P^yVb'y ) (106) 

where the baryon-to-photon ratio is given by 

^ ^ Pb + Pb _ 3 Pb 
P P'y 4 p-y 


( 107 ) 
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In the tight-coupling limit, = V;,. 

Using the coupled energy and momentum conservation equations we can write 


4'^7 


[(1 + R)Vh^]' 


-^V-Vb^-ch', 

-- (1 + . 


(108) 

(109) 


Neglecting the time-variation of 'h and R (i.e., 'h' ~ i?' ~ 0, which is a reasonable 
approximation) one gets the following oscillator equation: 


+ (1 + R)^ 


1 


3(1 -|- R) 


V" 


-5^ + (1 + 


0 . 


( 110 ) 


Comparing with the earlier oscillator equation neglecting baryons, (97), we see that 
the adiabatic speed of sound (the term multiplying the Laplacian is the square of the 
sound speed), and now is reduced by the presence of baryons by a factor (1 -|- i?). The 
solution for the above equation is: 


+ (1 + R)^ 


-d-y + (1 + R)'^ 


cos{ks) , 


( 111 ) 


J 0 


where now the sound horizon is 


s = / Csdr] = 


drj 


3(1 -|- R) 


( 112 ) 


So the matter-era solution for k ^ keq, (100), corrected by the presence of baryons is 
now: 


-I-PS [3i? — (1-I-3ii) cos(/cs*)]-C 7 • (113) 

The oscillations are now asymmetric: compression and rarefaction are no longer 
symmetric because baryons tend to collapse under their own gravity, more than they 
become more rarefied. 

In figure we present how the baryon loading modifies the shape of the CMB 
temperature power spectrum. We vary the baryon density but keep all the other 
parameters fixed. Because of Friedmann constraint ([^ 


1 - “t“ kliQ H/j; , 


(114) 
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Figure 3. CMB temperature power spectra, for different values of Vlbh?, showing the effect of 
baryon loading. We used the CAMB code m for the fiducial ACDM model given with the 
following parameters fixed = 0.022, VLch"^ = 0.12, = 0.7, r = 0.1 and flk = 0. 


we must therefore let h vary while keeping 0;, fixed. Increasing the baryonic density 
increases the height of the first peak and lowers the second peak. This is due to the 
fact that, the first peak is a maximum of compression, it is enhanced by a heavier 
baryon load, while the second peak is a rarefaction peak. 

In fig. I^we present how the dark matter content modifies the shape of the CMB 
power spectrum. We use the same strategy as in fig. but now we let flch? vary. 
As one can see, increasing this quantity causes the whole peak structure to decrease 
in amplitude. Since is fixed, the matter-radiation equivalence epoch takes place 
at earlier times when increasing flch?. This implies weaker radiation driving and a 
smaller amplitude of oscillations. 
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Figure 4. CMB temperature power spectra, for different values of showing the effect of 
radiation driving. 


There is also another very important feature characterising CMB anisotropies on 
small scales. On scales comparable to the photon diffusion scale the tight-coupling 
approximation breaks down. The diffusion comoving scale has the following expression 
in terms of cosmological parameters |10] 


Ad 


64.5 Mpc 


V 0.14 


- 0.278 


nbh‘^\ ' 
0.024 ) 


(115) 


Photon diffusion translates into a damping of the oscillations, see figure and a decay 
in the correlation, i.e., the angular power spectrum, Ci. This is easily understood. On 
very small scales, below Ad, cold and hot photons mix thereby averaging to zero the 
correlation. 

Thus far we have assume that all photons last scattered at the time of recom- 
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bination, but this is an approximation. The absence of a Gunn-Petersen trough (no 
absorption by neutral hydrogen in quasar spectra) suggests that neutral gas has been 
reionised by a redshift z > 6. Reionisation leads to an optically thin “smog” between 
us and recombination with an optical depth r ~ 0.1, i.e., 10 percent of the photons are 
scattered again. This reionisation suppresses small-scale anisotropies by rescattering 
the photons which also tends to average out the temperature anisotropies. The effect 
of varying the reionisation optical depth is shown in figure 



Figure 5. Effects of reionisation on the temperature (solid lines), EE (dotted lines) and BB 
(dashed lines) power spectra. 
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D. Parameter constraints from peak structure 


The dependence of the acoustic peak structure on a variety of different physical 
parameters enables cosmologists to determine these cosmological parameters with 
unprecedented accuracy from detailed measurements of the CMB sky. This is often 
referred to as the era of “precision cosmology”. 

From Planck, we have a precise constraint on the angular scale of the first acoustic 
peak [23]: 


(9* = (1.04148 ± 0.00066) x 10 


-2 


(116) 


In flat space this angle is given by 0* = i.e., the ratio between the sound 

horizon at recombination 


r* 2 

s* = / Csdrj oc ilmh , 

Jo 

and the present horizon distance to recombination 

dz 

/o Wz 

Thus in flat space we find 


D* = fjo — V* = / WWW oc h~^ 
Jo 


e^ = oc. nmh^ = 0.0959 ± 0.0006 . 


(117) 


(118) 


(119) 


More generally, in a curved space with curvature radius R, we have 0* = s^^/Da, 
see figure]^ The angular diameter distance. Da, can be written for <C 72 as 

n.mDV 


Da = 72 sin ( — 
from which we obtain the bound 




( 120 ) 


life = -0.0042lO;0^3 _ 


( 121 ) 


As yet there is no evidence for spatial curvature. 

We can derive a lot more information from the peak structure. For example, as 
discussed earlier, the second peak height relative to the first peak is related to the 
baryon loading. As shown in figure this suppresses the even (rarefaction peaks) 
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Figure 6. Effects of varying the cosmological constant (solid lines showing = 1) and 

spatial curvature (dashed lines showing flm + 11^ = 1). 

peaks with respect to the odd numbered peaks (compression peaks). From the relative 
height of the second peak, the following constraint on the the baryon content is found 

m- 


= 0.02207 ± 0.00033 . (122) 

The decay of the gravitational potential in the radiation era enhances the third and 
higher peaks. The less matter there is, the higher the peaks are enhanced, as shown 
in figure This enables us to put a constraint on the matter density today |24j : 

= 0.1423 ± 0.0029 . (123) 


Note that combining the two bounds (119) and (123) yields a bound on the Hubble 
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constant directly from the CMB 

Ho = 67.3 ± 1.2 km s“^ 


Mpc 


-1 


(124) 


The diffusion length (115) and hence the damping tail is then fixed once 
and are specified in the basic ACDM cosmology. Reionisation also suppresses 

anisotropies at all small angles, above I = 20. This is approximately degenerate with 


the primordial spectral tilt (85), but this degeneracy can be broken by polarisation 


(see next section). Planck data combined with WMAP polarisation data requires a 
spectral tilt 


n = 0.9603 ± 0.0073 . 


(125) 


VI. POLARISATION 

Thomson scattering is due to the motion of charged particles responding to 
an incident electromagnetic wave. The outgoing radiation emitted by an electron 
responding to a single incident wave is polarised in the direction of motion of the 
electron. However for the CMB photons emitted in a given direction at last scattering 
to have a net polarisation requires the electrons on the last scattering surface to see 
an anisotropic distribution of incoming photons. We have seen that tight-coupling 
between electrons and photons leads to an approximately isotropic distribution of 
photons before recombination, therefore CMB radiation from last-scattering is only 
weakly polarised, due to a small quadrupole moment caused by small, but finite photon 
diffusion before last scattering. However about 10% of CMB photons are re-scattered 
long after recombination due to reionisation, when the photon field is anisotropic, and 
hence reionisation provides a source of polarisation on large angular scales. 

In general we define the polarization tensor in terms of the electric field, E jinj 

, / 0-k(5 U \ 

Pii oc (EEt) oc (126) 

-u e-Q j 

where Q and U are the two possible states of linear polarization and the angle brackets 
here denote the time average. Circular polarization is neglected here. 
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Since the polarisation state is invariant under a 90° rotation, it corresponds to a 
spin-2 field. Like the temperature anisotropy, we can decompose the polarisation in 
any direction into harmonic functions, E and B modes, across the whole sky, 

(5(h) ± iU{n) = ^ {Eim ± iBim) ± 2 E£m(h), (127) 

£,m 

where ± 2 l£m(h) are spin-2 spherical harmonics. 

Scalar perturbations are longitudinal wave-modes where the inhomogeneities which 
give rise to polarised radiation are in the same direction as the wave propagates. This 
symmetry ensures that only E-mode polarisation is generated by scalar perturbations 
at linear order, corresponding to polarisation parallel or perpendicular to density 
gradients. Eigure shows the angular power spectrum for E-modes alongside the 
temperature power spectrum, and their cross-correlation. It shows the E-mode 
polarisation “bump” at large angular scales {i < 20) which can be used as a sensitive 
measure of reionisation and optical depth r. 

B-mode polarisation can only be generated at first order by transverse waves, i.e., 
vector or tensor perturbations. Initial vector perturbations decay at linear order 
in an expanding cosmology, and are completely absent in inflation driven by scalar 
fields. However tensor perturbations correspond to free perturbations of the metric, 
i.e., gravitational waves. 

Consider a spatial metric perturbation: 

(2^)3/2 ^ [hkeij(k) + hkeijik)] . (128) 

Tensor perturbations must be transverse and traceless, i.e., 

= 0 , Ecij = 0 , = 0 , = 0 . (129) 

This leaves only two independent degrees of freedom, which are the two polarisation 
modes characterising a gravitational wave. 

The Einstein equations give the following evolution equation for the gravitational 
wave amplitude 

P 

hk + 3Hhk H—= 0 , 


(130) 
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Figure 7. Effects of varying tensor-to-scalar ratio, r. The solid lines (top) correspond to the 
total intensity (temperature) angular power spectrum. The dashed lines correspond to the 
E-mode power spectrum and the dotted lines correspond to the temperature-E-mode cross¬ 
correlation. These are dominated by scalar modes. The dot-dashed lines show the B-mode 
angular power spectrum for different values of the primordial tensor-to-scalar ratio, r. 


i.e., a wave equation for a massless field, including a damping term due to the expansion 
of the Universe. Quantum vacuum fluctuations in massless fields during slow-roll 
inflation generate an almost scale-invariant spectrum of primordial perturbations on 
super-Hubble scales k < aH. Therefore, we can predict a primordial power spectrum 
for gravitational waves from inflation: 

327r ^2 


A /jr 

'p^^k) = 2-^{hh^) 


^Mj,i 2'k 


(131) 


k=aH 


{2tt)^ 

where the average is the vacuum expectation value, classically promoted to a variance. 
The power spectrum is thus directly determined by the Hubble parameter, H, during 
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inflation. 

It is customary to introduce the tensor-to-scalar ratio: 

_ T^t 

n ■ 

From Planck constraints on the temperature power spectrum we obtain [23j, r < 0.11, 
and the energy scale at the end of inflation is constrained to be y < 2 x 10^® GeV. 
Nonetheless there is considerable interest in searching for B-mode polarisation in the 
CMB, either from future satellite experiments or dedicated ground-based experiments, 
as a direct signal of primordial gravitational waves. If an almost scale invariant, 
Gaussian distribution of primordial gravitational waves were discovered then these 
would surely be strong evidence for inflation, and our first evidence for the quantum 
nature of gravity. 



VII. THE NEXT ERONTIER IN CMB THEORY 


The observed GMB sky is remarkably uniform with temperature variations less 
than one part in 10^. These lectures have reviewed the analysis of the homogeneous 
“background” CMB sky and anisotropies in the CMB temperature and polarisation, 
modelled using first-order perturbation theory. As the precision of CMB experiments 
improves the next challenge in CMB theory may be to accurately model non-linearity 
in the CMB anisotropies, both at last scattering and along the subsequent line-of-sight. 

Non-linear interactions lead to departures from Gaussianity, evident in the CMB 
bispectrum, i.e., a non-zero correlation between different spherical harmonics. Prim¬ 
ordial non-Gaussianity is often described in terms of the dimensionless nonlinearity 
parameter [38] 


Inl 


Bc_{ki,k2,k^) 

Pc{ki)Pc{k2) + Pc{k2)Pdk3) + Pc(fci)^c(^3) ’ 


(133) 


i.e., the primordial bispectrum, ^ 2 ,/cs), relative to the square of the power 

spectrum. This f^L ~ 1 corresponds to a primordial bispectrum B(^{ki, k 2 , k^) ~ 
10“^®. It is related to the three-point correlation function, which is identically 
vanishing in the Gaussian case, along with all the odd-order correlation functions. 
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In general Jnl defined in this way is a scale- and shape-dependent function of the 
three wave numbers, but for local-type non-Gaussianity Jnl is a constant parameter 

m- 

Non-linear interactions also lead to important conceptual differences from the 
simple assumptions valid in linear theory. Small but non-zero vector and tensor 
perturbations, hence B-mode polarisation, are generated at second-order from first- 
order scalar perturbations [IQI E]. Also, second order anisotropies in the the 
photon distribution can no longer be described simply by a black-body spectrum with 
anisotropic temperature. Second-order effects lead to anisotropic spectral distortions. 
For example, the angular power spectrum of the Compton y-distortion could provide 
a powerful tracer of reionisation |42j . 

Weak leasing along the line of sight is already an important nonlinear effect seen in 
current data, and needs to be taken into account in parameter estimates using Planck 
data. It is caused by many small-angle deflections by non-linear structures along the 
line of sight. This redistributes power in the small-scale angular power spectrum, 
smoothing out peaks at high i in the power spectrum. This anisotropic lensing of the 
small scale power in the CMB has been used by the Planck team to create a map of 
the lensing potential and hence a measure of the integrated mass distribution along 
the line of sight |43j . 

CMB weak-lensing also provides a non-zero contribution the angular bispectrum. 
Gravitational lensing leads to a second order anisotropy which is correlated with 
the integrated Sachs-Wolfe effect along the line of sight [44H46| . This contributes 
a significant bias to estimates of the non-linearity parameter, and /nl ~ 7 has been 
seen in the Planck analysis. After subtracting this effect, Planck measurements remain 
consistent with vanishing primordial non-Gaussianity, /tvl = 2.7 ± 5.8 HZ!- 

Most current bounds on primordial non-Gaussianity are based upon theoretical 
templates based on non-linear modelling of inflationary (or alternative) models, as 
non-Gaussian initial conditions for standard, linear Boltzmann codes such as CAMB 
133 or CLASS [IS]. However the process of decoupling and the Sachs-Wolfe effect on 
temperature (and spectrum) anisotropies is in reality a non-linear process. As bounds 
on primordial non-Gaussianity become tighter we also need templates for the intrinsic 
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non-Gaussianity expected from non-linear physics at recombination. There are now 
second-order general relativistic Boltzmann codes which have been developed |49H51j 
building on pioneering early work [52]. Intrinsic non-Gaussianity at last-scattering 
provides a small bias, /nl — 1; which remains below the observational uncertainty 
of current experiments. However we are now in the position to be able to build a 
template for the intrinsic non-linearity at recombination which could be a target for 
future all-sky (hence space-based) missions, as shown in figure This would make 
novel tests of physical process at last-scattering, e.g., gravitational wave production 
from density waves at second order. 



Figure 8. Signal-to-noise ratio of the intrinsic bispectrum signal as a function of maximum 
multipole for Planck and proposed future satellite missions using temperature only (dashed 
lines) or temperature and polarisation (solid lines). Figure reproduced from |53j . Copyright 
(2014) by The American Physical Society. 
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VIII. OUTLOOK 

The standard ACDM cosmology provides a remarkably successful base model, able 
to explain many detailed features of the CMB revealed over the past fifty years by 
a series of ground-, balloon- and space-based experiments. In particular the Sachs- 
Wolfe plateau at large angular scales, the series of acoustic peaks in the angular power 
spectrum, and the damping tail at small angular scales can be described by this model 
with just six cosmological parameters: the Hubble constant, the baryon and matter 
densities, the reionisation optical depth, and the amplitude and tilt of primordial 
perturbations. These six parameters are increasingly tightly constrained in the new era 
of precision cosmology, and bounds are set to become ever tighter, especially through 
new combinations with other data sets, such as high redshift galaxy surveys and HI 
(neutral hydrogen) survey data. The framework already successfully accommodates 
new observational discoveries such as the effect of weak lensing now seen in the CMB 
power spectrum and bispectrum. 

Nonetheless there is no reason to believe this is the final theory of cosmology. 
Even the simplest, single-field models of inflation in the very early universe make 
predictions for additional features in the primordial perturbations, including a 
spectrum of tensor (gravitational wave) perturbations and small but finite running 
of the scalar spectral index. Many inflation models make further predictions including 
primordial isocurvature perturbations and/or non-Gaussianity. Any of these additional 
parameters would radically change our views about the likely mechanisms generating 
primordial structure. There are many additional cosmological parameters possible, 
including additional particle species and/or interactions, but there is no clear evidence 
yet requiring any more than the six basic parameters. 

The present theoretical framework now being constrained by data was established 
in the 1970s and 80s well before the golden age of CMB experiments was begun 
by the COBE satellite results. Work now in progress will set new theoretical 
challenges for future experiments. Ground-based experiments are currently targetting 
CMB polarisation from weak-lensing and the elusive B-mode signature of primordial 
gravitational waves. Future space-based experiments are likely to focus on polarisation. 
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spectral distortions and/or non-Gaussianity. The CMB will remain a key testing 
ground for cosmological theory for many years to come. 
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